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Measurements of the hadronic final state at HERA are reviewed, which aim at the 
investigation of the parton dynamics of the proton at small Bjorken x. 



1 Introduction 

Amongst the most interesting issues of HERA physics is QCD in the newly 
accessible regime of small Bjorken x. The observed rise of the structure func- 
tion i<2 towards small x suggests a strong increase of the parton density in the 
proton, but what is its dynamical origin? Is BFKL |jj dynamics at work, or is 
conventional DGLAP || evolution sufficient? Complementary measurements 
of the hadronic final state provide more detailed information than the inclusive 
i*2 data to help uncover the underlying dynamics. 

The leading log DGLAP resummation corresponds to a strong ordering of 
the transverse momenta kx (w.r.t. the proton beam) in the parton cascade, 





Fi gure 1 '. a) Ladder diagram for parton evolution, b) The multiplicity of hard gluons 
with p T > 2 GeV vs. CMS n. The events are generated with the models CDM, MEPS 
and HERWIG in a "low x " kinematic bin with (x )=0.00037 and {Q 2 >=13.1 GeV 2 . The 
proton direction is to the left. 
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Qo < k T1 < — ^Ti < ---Q 2 (fig- a )- Since in the BFKL evolution that re- 
striction is absent ( "fcr diffusion" ) , a generic signal for deviations from DGLAP 
evolution is enhanced radiation from the ladder between the current and the 
remnant system, that is in the central rapidity region in the hadronic centre 
of mass system (CMS). Experimentally pseudorapidity rj = — lntan(0/2) is 
used, where 9 is the angle with respect to the virtual photon axis. The fol- 
lowing observables have so far been exploited and are being discussed in this 
paper: 1) Et flows: increased parton activity should result in an increased 
transverse energy flow 2) Charged particle pt spectra: high kx partons, 
disfavoured by the strong kx ordering in DGLAP, are signalled by measure- 
able high pt hadrons 3) Forward jets: high energy jets with prjet 2 ~ Q 2 
(kinematically bound to be measured in the forward calorimetric systems close 
to the remnant direction) tag events with BFKL evolution, because DGLAP 
evolution is not allowed |5|, |(| . 

Theoretical calculations for these observables exist and can be compared to 
the data, provided that hadronization effects can be controlled. Alternatively, 
predictions are derived from Monte Carlo models, which incorporate the QCD 
evolution in different approximations and utilize phenomenological models for 
the non-perturbative hadronization phase. The MEPS model (Matrix Element 
plus Parton Shower, program LEPTO @) and the HERWIG generator || are 
based upon leading log DGLAP parton showers, with strong kx ordering of the 
emitted partons. In the colour dipole model (CDM) ||, [H)| gluon emission is 
not subject to kx ordering [pjj . In that respect it mimicks the BFKL evolution, 
and leads to more abundant gluon radiation than in the other models (fig. |l]b). 



2 Energy flows 

The flow of transverse energy Et as a function of 77 provides a very simple, 
global characterization of the hadronic final state. Though their partonic ra- 
diation patterns are very different, all models in their present incarnations 
provide a reasonable description of the detailed Et flow data by HI []l2| (two 
out of 17 x — Q 2 bins are shown in fig. ^). 

From BFKL evolution a relatively large amount of Et at central rapidity is 
expected that increases with decreasing x, opposite to the DGLAP expectation 
. These behaviours are realized at the parton level for the models with the 
unordered and ordered emission scenarios (fig. ^). The average Et ((Et)) 
measured in —0.5 < r\ < 0.5 |l2| , [H| does increase with falling x, in qualitative 
agreement with the BFKL expectation (fig. |^). However, when comparing to 
the calculation hadronization effects have to be taken into account. 

The DGLAP based models have to employ large hadronization corrections 
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Figure 2: The E T flow vs. r) in the hadronic CMS for (x) ss 0.002 and (x) RJ 0.0004 with 
(Q 2 ) fa 14 GeV 2 fixed. The proton direction is to the left. The data Jl2| are compared to 
the models CDM (ARIADNE 4.08), MEPS (LEPTO 6.4), and HERWIG 5.8. 



to achieve the level of Et seen in the data (see fig. ||b) . For example, in LEPTO 
the new concept of soft colour interactions (SCI) has been introduced, resulting 
in a reconfiguration of the fragmenting strings, that may lead to an enhance- 
ment of Et (and also to rapidity gaps) (li). In the CDM, hadronization effects 
are much smaller than in LEPTO and HERWIG (fig. ||). In summary, though 
the data follow the trend expected from BFKL evolution and are consistent 
with a model with an unsuppressed radiation scenario, they are also consistent 
with DGLAP evolution, assuming large hadronization effects. 

3 Charged particle transverse momenta 

Not yet well understood hadronization effects precluded strong conclusions 
on the underlying parton dynamics from the Et flow measurements. Single 
particle transverse momentum (py ) spectra represent a more direct measure of 
the partonic activity [0 . Unsuppressed parton radiation should manifest itself 
in a hard tail of the pt distribution, whereas hadronization should produce 
typical spectra limited in pt ■ 

HI has measured the charged particle pt spectra as much central in the 
CMS (0.5 < rj < 1.5) as the tracking detector acceptance allowed J^]. At 
large x, all models agree with the data, but at small x, only the model with 
the unsuppressed radiation pattern (CDM) is able to describe the high pt tail 
seen in the data (fig.[|a). The shortfall of the models with suppressed gluon 
radiation indicates that at small x there is more high Ut parton radiation 
present than is produced by the models based upon leading log DGLAP parton 
showers. 
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Figure 3: The mean £?r (GeV) in —0.5 < rj < 0.5 as a function of x for Q 2 = 14 GeV 2 . 
The data |l| are shown together with the models CDM (ARIADNE 4.08), MEPS 
(LEPTO 6.4), and HERWIG 5.8 for hadrons and for partons, and with the BFKL cal- 
culation for partons The LEPTO result for hadrons without the features soft colour 
interaction (SCI) and the new sea quark treatment (SQT) is also shown. 



The pt spectra at small x have been calculated by folding a cross section 
aj to produce a parton j with experimentally known fragmentation functions 
Df,/j(z) to produce a hadron h with momentum fraction z from the parton j 
pa. Symbolically, er^ = tTj^D^/j. Monte Carlo assumptions for hadronization 
are thus avoided. When BFK1 evolution is invoked, the HI data are well 
described by the calculation (fig |^b) . The normalization for the BFKL part 
is obtained by requiring that the calculated parton cross section would match 
the measured HI forward hadron jets |i~7j (see next sect.). Neglecting BFKL 
evolution, the calculation falls significantly below the data. It would now be 
interesting to compare a complete fixed order calculation in next-to-leading 
order (NLO) to the data. 



4 Forward jets 

Forward jets are the classic signature for BFKL evolution |(|. One requires 
Xj ct = Ej e t/Ep, the ratio of jet and proton energy to be large in order to maxi- 
mize the BFKL evolution from Xj C t down to x. The jet transverse momentum 
^Tjethas to be close to Q, fcrjot ~ <5 to suppress the phase space for DGLAP 
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Figure 4: a) The p T spectra of charged particles from 0.5 < r\ < 1.5 (CMS) |15| . 
Displayed are two different kinematic bins at high and low x for (Q 2 ) ~ 14 GeV 2 . The 
models ARIADNE 4.8, LEPTO 6.4 and HERWIG 5.8 are overlayed. b) The result of a 
theoretical calculation with and without BFKL evolution [juj, compared to the HI data at 
low x for different choices of the factorization scale /i. 



evolution. In the presence of BFKL evolution, the forward jet cross section 
should grow faster with decreasing x than for DGLAP evolution. 

HI reconstructs forward jets with the cone algorithm and requires prjet > 
3.5 GeV p7| . The forward jet cross section (fig. ||a), corrected for detector 
effects, is well described by the CDM, and somewhat less well by the standard 
MEPS model with SCI (note that SCI was necessary in order to describe the 
E T flows). Without SCI, the DGLAP based model MEPS cannot describe the 
growth of the cross section towards small x. Comparing the hadron level jets 
to the parton level jets (fig. ||b), one finds relatively small hadronization correc- 
tions in the CDM and also in MEPS without SCI, but with SCI, hadronization 
effects become large! If one would assume the CDM hadronization corrections 
to be correct, the data would be far above both the NLO calculation |l^] and 
the parton level jets from the DGLAP based model MEPS. 

Presumably hadronization uncertainties become smaller for a higher PTjet 
cut-off. In the ZEUS analysis j^O), which relies also on the cone algorithm, 
PTjet > 5 GeV is required. The data in fig. |^c are corrected to the parton 
level, assuming hadronization corrections from the CDM. In contrast to a NLO 
jet calculation plj , the CDM describes the data well. 

It is clear that in order to draw firm conclusions from the forward jet 
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Figure 5: a) The HI forward jet cross section |17| (1994 data) vs. x , corrected for detector 
effects to the hadron level. Xj et > 0.035, 0.5 < PTjct 2 /Q 2 < 2 and PTjct > 3.5 GeV were 
required. Also shown are the hadron level predictions from CDM (ARIADNE 4.08) and 
MEPS (LEPTO 6.4, with and without SCI), b) The corresponding parton level forward jet 
cross section from a NLO calculation [[li|, and for CDM and MEPS. Also shown is a BFKL 
parton cross section calculation Jl9| (no jet algorithm), c) The ZEUS forward jet cross section 
pp[ (1994 data) vs. x , corrected to the parton level. Xj et > 0.035, 0.5 < PT)ct 2 /Q 2 < 4 
and PTjct > 5 GeV were required. The data are compared to a NLO jet calculation [ plj 
and to parton jets from CDM (ARIADNE 4.08). Also shown are parton cross sections (no 
jet algorithm) with ("BFKL") and without ("Born") BFKL evolution Jis| . The systematic 
errors do not include uncertainties due to hadronization. 



analyses, hadronization effects need to be better understood, ft is however 
interesting to note that the BFKL calculation |l9j is far above the Born level 
calculation (excluding BFKL evolution) , and describes qualitatively the rise of 
the measured forward jet cross sections at small x (fig. ^|). Direct comparison 
to the data however would be imprudent, because these calculations do not 
invoke a jet algorithm, and hadronization corrections are potentially large. 



5 Conclusion 

The transverse energy flow and the forward jet data are only compatible with 
a conventional DGLAP evolution scenario when large hadronization effects are 
assumed. The high px tail seen in the charged particle spectra cannot be ex- 
plained with hadronization effects. A theoretical calculation neglecting BFKL 
evolution, where hadronization is taken into account via known fragmentation 
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functions, falls far below the data. The pr spectra require more hard parton 
radiation than expected from conventional DGLAP evolution. BFKL effects 
offer a consistent explanation of the measured px spectra, forward jets and 
Et flows. Of course this does not exclude other explanations. For example, it 
has been suggested that contributions from resolved virtual photons could be 
responsible for enhanced parton activity [£2[ . In that case the hard scattering 
takes place between partons inside the resolved photon and partons from the 
proton and may thus happen not at the photon vertex but further down the 
ladder. 
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